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ABSTRACT 

Under nonlocal convection theory, convection extends without limit therefore no ap- 
parent boundary can be defined clearly as in the local theory. From the requirement of 
a similar structure for both local and non-local models having the same depth of con- 
vection zone, and taking into account the driving mechanism of turbulent convection, 
we argue that a proper definition of the boundary of a convective zone should be the 
place where the convective energy flux (i.e. the correlation of turbulent velocity and 
temperature) changes its sign. Therefore, it is convectively unstable region when the 
flux is positive, and it is convective overshooting zone when the flux becomes negative. 
The physical picture of the overshooting zone drawn by the usual non-local mixing- 
length theory is not correct. In fact, convection is already sub-adiabatic (V < Vad) 
far before reaching the unstable boundary; while in the overshooting zone below the 
convective zone, convection is sub-adiabatic and super-radiative (Vrad < V < Vad)- 
The transition between the adiabatic temperature gradient and the radiative one is 
continuous and smooth instead of a sudden switch. In the unstable zone the temper- 
ature gradient is approaching radiative rather than going to adiabatic. We would like 
to claim again that, the overshooting distance is different for different physical quanti- 
ties. In a overshooting zone at deep stellar interiors, the e-folding lengths of turbulent 
velocity and temperature are about 0.3Hp, whereas that of the velocity-temperature 
correlation is much shorter, being about 0.09Hp. The overshooting distance in the 
context of stellar evolution, measured by the extent of mixing of stellar matter, should 
be more extended. It is estimated as large as 0.25-1.7 Hp depending on the evolution- 
ary timescale. The larger the overshooting distance, the longer the timescales. This is 
due to the participation of extended overshooting tail in the mixing process. 
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1 INTRODUCTION 

As the classical treatment of convection, the local theory 
has been used in modelling stellar structure and evolution. 
In the calculation of massive star evolution, Schwarzschild 
& Harm (1958) discovered that the hydrogen rich radiative 
envelope just outside the helium rich convective core cannot 
be convectively stable, and that leaded to the paradox of 
so called semi-convection. To solve that problem, the idea of 
semi-convection was initiated, i.e. the region outside the con- 
vective core is in a state of semi-convection. Stellar matter 
in this region is nearly in neutral stability (V < V ad), there- 
fore convective energy transport due to this mild convection 
can be neglected, while the mixing of chemical compositions 
should be important, which makes a gradient of molecu- 
lar weight in this region (otherwise called semi-convection 
zone). There had been a great debate in the community for 



a long period since then on whether the Schwarzschild or 
Ledoux criteria should be applied for the neutral stability 
of convection, and whether the semi-convective zone should 
be very wide or rather narrow. Stothers (1970) commented 
on various establishments of semi-convection. Evolutionary 
scenarios for massive stars with or without semi-convection 
were also discussed (eg. Chiosi & Summa 1970). It has been 
realized later that the problem of semi-convection is in fact 
due to the non-locality of stellar convection. Therefore var- 
ious theories of non-local theory of stellar convection have 
been worked out (Spiegel 1963, Ulrich 1970, Xiong 1977, 
1981a, 1989, Kufuss 1986, Grossman et al 1993, Canuto 
1993, Canuto & Dubovikov 1998). Such non-local theories 
of stellar convection have then been applied in the studies 
of structures of the solar and stellar convective envelopes 
(Travis & Matsushima 1973, Unno et al. 1985, Xiong & 
Cheng 1992), stellar oscillations (Xiong 1981b, Xiong, Deng 
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& Cheng 1998, Xiong, Cheng & Deng 1998, Xiong & Deng 
2001, 2007) and stellar evolution (Xiong 1986). Generally 
speaking, non-local theory of convection makes the results 
better match observations than the local ones. However, the 
non-local theory of convection is rather complicated, which 
is much less straightforward to be understood, much more 
difficult to be applied and much more computing power de- 
manding than the phenomenological local (Bohm-Vitense 
1958) and non-local mixing length theories (Maeder 1975, 
Bressan et al. 1981). For these reasons, it becomes a gen- 
eral practice to use the phenomenological local or non-local 
mixing length treatment for stellar convection in nowadays 
stellar evolution models. The non-local mixing of chemical 
compositions during the evolution of stars is dealt with by 
attaching a parametric overshooting zone outside the con- 
vectively unstable region. The parametric distance of con- 
vective overshooting has a great impact on the properties of 
stellar evolution. The goal of present work is to discuss the 
calibration of the overshooting distance. The physical defi- 
nition of the boundary of the convective zone is discussed in 
the next section. In section |3l calibrations of the overshoot- 
ing distance by numerical simulations of non-local convec- 
tion and depletion of solar Lithium abundance are presented. 
A summary and discussions are given in the last section. 



2 HOW TO DEFINE THE BOUNDARIES OF A 
CONVECTIVE ZONE? 

Normally in a local theory of convection, the boundary of 
convective zone is given by the so called Schwarzschild's cri- 
terion, 

V = V,d, (1) 

which is derived by analysis of the local convective sta- 
bility. However, if viewed from hydrodynamics strictly, all 
hydrodynamic phenomena including convective motions in 
stars are non-local, therefore there should be no well defined 
boundary for convective motion in an extended medium. In 
a certain sense, forcing a definition of boundary for a convec- 
tive zone is always an artifact. Defining a boundary for stel- 
lar convection is needed in practice of stellar evolution cal- 
culations, but this cannot be done arbitrarily, instead some 
objective standards should be respected. These standards 
should at least include the following: 

(i) As a matter of fact, most of calculations for stellar 
structure and evolution still use local theory of convection. 
Therefore, the definition of boundary given by a non-local 
convection theory should be kept as close as possible to that 
by a local theory. In other words, the local and non-local 
convection models with the same depth of convective zone 
should be made to have structures as similar as possbile. 

(ii) the definition of the boundary should be physically 
pounced in any case, i.e. the unstable convective zone should 
be the driving (excitation) region of convective motion, and 
the adjacent overshooting zone should be the dissipation 
region of convective motion. 

Stellar convection happens due to some internal insta- 
bility of the thermal structure in gravitationally stratified 
fiuid, therefore the study of the resulted convective motions 



should be based on the dynamical equations of fluid. A com- 
plete dynamic equations of time-dependent non-local con- 
vection theory can be found in our previous work (Xiong 
1981, 1989). For the sake of clarity, and to make it easier to 
read this paper, the derivation of the dynamic equations of 
turbulent convection in steady fiuid is presented here. The 
conservations of momentum and energy of fluid dynamics 
can be expressed as. 



d IpU J , i k , ik 7->\ 

+ Vk [pu u + g P) 



dt 



+ Vfc [pu H) 



dt 



dP 
'dt 



u^VkP + ^kFr = pen + (T'''(u)VfeMi 



(2) 



(3) 



where p, T and P are the regular labels for density, 
temperature and pressure (including radiative pressure) of 
gas; H and axe enthalpy and nuclear energy generation 
rate per unit mass; is the ith component of fiuid motion 
vector; cr'*(ij) the viscous stress tensor; is ith component 
of the radiative fiux vector. The implicit summation rule 
of tensor calculations is used, i.e. a pair of sub- and super- 
script index mains summation from 1 to 3. When convection 
happens, any physical quantity can be written as the sum 
of averaged and turbulent fiuctuated components as, 
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Putting the expressions of all quantities in the form of 
eq. @ into eqs. ((2}-([3]), and averaging the whole equations, 
the dynamic equations for the mean fiow can be derived as. 
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Subtracting the corresponding equations of the mean 
motion eqs. ©-(O from eqs. ©-(IS}, and considering the 
static state of the flow, i.e. 

_ dP dp dH 

the dynamic equations for the fluctuation quantities can 
be derived as. 
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By using a certain thermodynamic relations, and fol- 
lowing some deductions and simplifications, eq. ((9} can be 
written as, 
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where w' is the density weighted fluctuation of turbulent 
velocity, 
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Starting from eq. (|8]) and eq. (|10p . we can have 
the following dynamic equations for the auto- and cross- 
correlations of turbulent velocity and temperature fluctua- 
tions: 
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where , Z and V axe respectively the auto and cross 
correlations of turbulent velocity w' and the relative tem- 
perature fluctuation T' /T, defined as the following, 



(14) 



x^ = w'^w'JS, 
Z^{T'/f)\ 



V = w'rT'/T, 



(15) 
(16) 
(17) 



while V ad is the adiabatic temperature gradient, V = 
dlnT /din P is the temperature gradient, Xc is a variable 
related with the effect of thermal conductivity: 
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Pe — x/xc is the effective Peclet number of turbulent 
convection. B = — {din p/dlnT) p the expansion coefficient 
of gas. The detailed derivation of the dynamic equations 
of correlations can be found in our previous work (Xiong 
1978, 1981, 1989). Eqs. ([T2l),([T3l) and ^ are the dynamic 
equations of turbulent convection in steady fluid, which have 
very clear physical meanings. Eq. (|12p . for instance, is for 
the conservation of turbulent kinetic energy. The left hand 
side is the rate of variations of turbulent energy per unit 
volume, which is equal to the sum of the 3 terms on the 
right hand side: the first term is the work done by buoyant 
force: 



(19) 



while the expression in the bracket of the second term 
represents the effect of non-locality, which is the fiux of tur- 
bulent kinetic energy Lt, 



Lt = 4nr pu'r'w'rW'^/2. 
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Therefore the second term on the right hand side of 
eq. (|12p is the net gaining rate of turbulent kinetic energy 
per unit volume —pdLr/dMr- The third term on the right 
hand side of eq. (|12|) is for viscous dissipations, i.e. the dis- 
sipation rate pel by converting turbulent kinetic energy into 
thermal energy due to viscosity. Ci is a convection parameter 
related to the viscous dissipation of turbulent convection. 
le = ciHpr / Ro is the linear size of the energy-containing 
eddies (Xiong 1978, 1981a, 1989). Both Eqs. ^ and dl 
have similar physical means: the left hand sides are the vari- 
ation rates of correlations Z (or V) respectively. The first 
(and the second for eq. [TJ) term on the right hand side 
of eq. (|13|) [ or eq. (|14|) ] is the rate of increase of Z (or V) 
due to super-adiabatic temperature gradient (and buoyant 
force for eq. dH); the second (the third for eq. [14] ) term 
represents the net increase rate of Z (or V) due to non-local 
convective energy fiux; while the last term describes the tur- 
bulent dissipation due to viscosity and thermal conductivity. 
For static convection, the time derivatives terms on the left 
hand side of eqs. (|12|| - (|14p all vanish: 

dx^ _dZ _dV 
~dt ~ ^ ~ ~dt 

Hence, eq. (fT2)) can be rewritten as, 
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It is clear from eq. ([22} that, for static convection, the 
net gaining rate of turbulent kinetic energy {-pdLr/dMr) 
and buoyant force work (Wbuo) will be balanced by turbulent 
dissipation pel (the sum of all the three terms vanishes). 
There is no need to go through the similar physical meanings 
of eqs. p3p and (|14p . When neglecting all the third order 
correlations representing the transportation effect of non- 
local turbulent convection, eqs. (I12p - (|14|l can be converted 
to, 



BV ~1.56-^x^ = 0, 

ClP 



(V - V,d) V - — {x + x,)Z ^0, 

Cl 

P 78 

(V - V„d) + B-Z - — {3x + Xc) V = 0, 

p Cl 



(23) 
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It is rather easy to solve eqs. ((23)) - (l25ll for x^ Z and V, 




Comparing eqs. (|26|) - (|28| ) and the equations of the local 
mixing length theory (Bhom-Vitense 1958), and substituting 
Ci and the effective Peclet number x/xc in our equations by 
the mixing length parameter a and 7 respectively, it is clear 
that the two expressions are the same. The stability condi- 
tion for convection is the Schwarzschild criterion (eq. [1]) in 
a chemically homogeneous medium. When there is a molec- 
ular weight gradient, the neutral stability condition should 
be Ledoux criterion (Xiong 1981a); 

V = V„d + Vp, (29) 

Hence, viewed from hydrodynamics, the local mixing 
length theory is only a special simplified case of our sta- 
tistical theory of correlations for turbulent convection. The 
third order correlation terms in eqs. (|12p - (|14p represent non- 
local effect of turbulent convection, neglecting which makes 
the equations becoming the local expressions eqs. (|23|l - (|25|) 
or their explicit form eqs. (|26p - (|28|l . In this case, a con- 
vective zone will have a clearly defined boundary given by 
Schwarzschild (or Ledoux) criterion: V > Vad is convec- 
tively unstable, while V < Vad stable (radiative). As we 
shall show later, within the convectively unstable [V > 0) 
zone far away from the boundary, the third order correlation 
terms in eqs. (|12|| - (|14p can be safely neglected compared 
with other terms. This means that the local expression of 
convection is a fairly good first approximation at the deep 
interior of an unstable zone. This is exactly the reason why 
the mixing length theory is still widely applied in the cal- 
culations of stellar structures. However, when studying the 
entirety of a convective zone, especially near the boundary 
of convective zone and in the overshooting region, the third 
correlation terms cannot be neglected. Instead, they are the 
true reasons for the existence of convective overshooting. In 
the non-local convection theory, it is clear from eqs. (|12|| - 
(|14p that the turbulent velocity and temperature fluctua- 
tions are different from zero everywhere, 

x > 0; Z >Q. (30) 

These conditions make it difficult and uncertain to de- 
fine a boundary of a convective zone. Schwarzschild criterion 
is overwhelmingly used in the community to fix the bound- 
ary. They think that the temperature gradient is very near 
and slightly higher than the adiabatic temperature gradient 
in the unstable region; and is also near but slightly lower 
than the adiabatic one within the overshooting zone. We 
are going to show that such a picture for convective over- 
shooting is not correct following the dynamical theory of 
turbulent convection. The cause of such a mistake is that 
there is a implicit hypothesis that has been applied in the 
phenomenological mixing length theory: turbulent velocity 
is fully correlated (either positively or inversely) with tem- 
perature (see eg. Xiong & Cheng 1992, Petrovay & Marik 



1995). In fact, when convection is very effective (the effective 
Peclet number x/xc S> 1), the correlation between turbu- 
lent velocity and temperature fluctuations near the bound- 
ary and in the overshooting zone decreases very quickly and 
vanishes eventually (Xiong & Cheng 1992). 

Eqs. (|12p -- (|14|l are derived under rather general condi- 
tions among which are two important assumptions as the 
following: 

(i) convection is subsonic, the relative fluctuations of tem- 
perature and density are both far less than unity: 

\p'/p\ « 1; \T'/f\ « 1, (31) 

(ii) inelastic approximation which actually filters out all 
acoustic waves not important for energy transfer in subsonic 
convection. 

It is well known that the dynamic equations of turbu- 
lent correlations have no closure due to the nonlinearity of 
hydrodynamics. That means: the third order correlations 
must be present in the dynamic equations of the second or- 
der correlations; while the fourth order ones bound to turn 
up in the equations of the third order ones, and so forth. 
Some hypothesis must be used in order to make a closure 
for the dynamic equations of the correlations. Obviously, the 
closure cannot be unique. Quite a few methods have been 
adopted so far (Xiong 1981a, 1989a, Canuto 1993, Gross- 
man et al. 1993, Canuto & Dubvikov 1998). In our opinion, 
a good closure should meet the following conditions: 

(i) the solutions of the resulted equations must be phys- 
ically sensible. For instance, the standard quasi-normal ap- 
proximation seems to be better in terms of mathematics, 
however it gives solutions like x^ < Q ox Z < Q which 
are physically non-sense (Grossman 1996). Such a seemingly 
reasonable assumption, if not modified somehow, cannot be 
used for the closure of the dynamical equations of the third 
order correlations; 

(ii) the solutions presented should not be in contradic- 
tion with observations. For instance, it should reproduce the 
main observational properties of solar granular velocity field, 
it should be able to explain the pulsation instabilities of low 
temperature stars having extended convective envelope, and 
it should be able to model the observed Lithium abundance 
patterns in the atmospheres of the Sun and solar type stars, 
and so on; 

(iii) the solutions provided should be comparable to that 
of direct hydrodynamical simulations. 

Our non-local theory of convection (Xiong 1981a, 
1989a) have been tested against the above standards, quite 
satisfactory results have been reached, therefore we have 
good reason to believe that it has nicely expressed the dy- 
namic behaviors of stellar turbulent convection. 

The solid line in fig. [1] shows the fractional convective 
flux Lc/L versus depth logP for a model of the solar con- 
vective zone calculated with our non-local convection theory, 
and the dashed line is that of a local convection model hav- 
ing the same depth of convective zone. It follows from the 
flgure that there is almost no difference between the local 
and non-local model, except some sizable deviations near 
the boundary of the convective zone; Fig. [5] depicts the rel- 
ative squared sound speed and density difference between 
the local and non-local solar convection zone models with 
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Figure 1. The super-adiabatic temperature gradient V — Vadi 
turbulent velocity-temperature correlation RvT, and the frac- 
tional convective flux Lc/L versus the depth (logP) for a non- 
local convection model of the Sun. The dashed line is the frac- 
tional convective flux L^/L for a local model with the same depth 
of convective zone. 

the same depth of convective zone versus depth. The rela- 
tive diflerence between the two models is mostly below 1% 
excluding the solar surface region. In this case, the bound- 
ary of the (non-local) convective zone is set at where the 
turbulent velocity-temperature correlation vanishes, 

V ^0; (32) 

Passing through the boundary, V changes its sign: 
within the convective zone: 

V >0, (33) 

and in the overshooting zone: 

V <0. (34) 

It is clear from Figs. [T] and [5] that if the boundary is 
defined as such, the structures of the local and non-local 
convection models with the same depth of convective zone 
should be similar. It can be understood by the fact that, 
within the stellar interior, the turbulent kinetic energy fiux 
(Lt) is generally much less than that of thermal convec- 
tion (Lc), and the turbulent pressure {Pt = px^) is much 
less than that of gas (Pg). It can be easily shown that 
Pt/Pg ~ x^/Cs = Ma^, where Cs is the local sound speed. 
Ma is the Mach number of turbulence; Lt / Lc is of the same 
order of magnitude as that of Pt/Pg. Except at the top of 
convective zone, we have Ma <^ 1. It follows from fig.|3]that, 
for the Sun, Lt/L < 1%, therefore the thermal convection 
Lc dominates the pressure-temperature (P-T) structure. It 
is then clear that, when defining the boundary of convective 
zone by V changing its sign, the structures of the local and 
non-local models having the same depth of convective zone 
will be similar. 
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Figure 2. The relative differences in the squared sound speed and 
density between the non-local and the local convection models 
with the same depth of the convective zone versus the fractional 
radius. 



Fig. [T] clearly demonstrates that convective motions 
near both upper and lower boundaries of the convective 
zone are very different. This is due to the fact that, in the 
atmosphere, the density is very low and Pe = x/xc < 1, 
therefore convective energy transfer is inefficient. As a re- 
sult, there exists a thin super-adiabatic layer atop of the 
convective zone. Passing through the upper boundary, the 
turbulent velocity-temperature correlation Rvt = V/xZ^^^ 
drops quickly from ~ 1 to — 1. This theoretical predic- 
tion agrees the observations of solar granular velocity field 
(Leighton et al.l962; Salucci et al. 1994) and the results of 
hydrodynamic simulations (Kupka 2003). Contrary to the 
situations in the solar atmosphere, convection is highly ef- 
ficient in terms of energy transfer (Pe ^ 1) in the deep 
interiors of the Sun. Towards the lower boundary of the con- 
vective zone, the turbulent velocity-temperature correlation 
Rvt decreases abruptly and approaches zero (|iivT| ^ 1). 
What makes it so different at the two boundaries is the dis- 
tinct the effective Peclet number. 

Fig. [5^ shows the work done by buoyant force pWbuo, 
the net gain (> 0) or loss (< 0) of turbulent kinetic energy 
due to non-localism of turbulence —pdLr/dMr and turbu- 
lent viscous dissipation rate pe as functions of depth in the 
solar convective zone. As shown in fig. [3^, contribution due 
to non-local convection is much less than the other two quan- 
tities in the convective zone, except in the narrow regions 
near the boundary of convective zone and in the overshoot- 
ing zones. The energy balance in the balk of convective zone 
is due to the interplay of the work done by buoyant force 
Wbuo and the viscous dissipation pe. Turbulence retrieves en- 
ergy from buoyant force, while at the same time it dissipates 
energy due to viscosity. The former factor, as the source, 
originates primarily from large eddies; while the later one 
is happening in the viscous dissipation range of the high- 
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Figure 3. a), the work done by buoyant force, viscous dissipation of turbulence pe, fractional turbulent kinetic energy flux Lt/L and 
the net gaining rate of kinetic energy due to turbulent diffusion —pdLt/dMr versus depth (log P) for a non-local solar model, b). and 
c). are the expanded plots near the upper and lower boundaries of convective zone indicated by the vertical dashed lines. 



est end of the turbulent spectrum: turbulence gains energy 
from buoyant force, which is then cascaded from lowest to 
higher and higher wave numbers of turbulent spectrum, and 
is eventually converted into thermal energy due to molecular 
viscosity. It follows from fig. [3^ that the work done by buoy- 
ant force pWbuo and the viscous dissipation amount nearly 
the same but have opposite sign within the convectively un- 
stable region. 

Figs.|3}3 and|3}; is the enlargements of fig.|3K around the 
two boundaries. Learnt from the two plots, attention should 
be paid to the following two points; 

(i) In the overshooting zone, both e and Wbuo are neg- 
ative, while dLr/dMr > 0. Therefore it is the non-local 
convective diffusion that drives overshooting, without which 
there would be no overshooting; 

(ii) Before reaching the boundary from the unstable side, 
the super-adiabatic temperature gradient has already be- 
come negative, which is more prominent near the lower 
boundary of the convective zone (as in fig.|3|. The boundary 
of convective zone {V = 0) is located at log P = 13.74 (the 
vertical dashed line), but at log P « 13.34 still in the unsta- 
ble zone, convection is already sub-adiabatic (V — Vad < 0, 
see also fig. [1}. This is distinctly different from the phe- 
nomenological non-local mixing length theories. Although it 
is resulted from our special non-local mixing length theory, 
such properties of convection ought to be general. This can 
be proved by eq. (|14p of the general dynamic equations of 
turbulent convection: The first two terms on the right hand 
side of eq. (|14p are both positive (a:, Z > 0). As discussed 
already, the third order correlations representing the non- 
locality of convection can be neglected compared with the 
second order terms in the deep interior of convective zone. 
Therefore, when V becomes zero (the fourth term on the 
right hand side), the super-adiabatic temperature gradient 
must be negative (V — Vads < 0) to make the equation math- 
ematically right. This proves that V — Vad must turn neg- 
ative before V does approaching the boundary. This nature 
of non-local convection does not depend on the kind of non- 
local convection theory used, which is universally shared by 
all theories of non-local convection following hydrodynam- 
ics. 



3 HOW EXTENDED IS CONVECTIVE 
OVERSHOOTING 

Due to the complexity of non-local convection theory, almost 
all the modellings of stellar structure and evolution are still 
using the local convection theory. Convective overshooting 
is defined as the penetration of convective motion through 
the classical boundary of convectively unstable zone into 
the adjacent stable region. The extent of overshooting is not 
the same for different physical quantities following our dy- 
namic theory of convection, and this leads to some troubles 
in understanding and estimation of overshooting. Followed 
by the great success of helioseismology, people are expecting 
to draw a firm conclusion to the long debated overshooting 
distance at the bottom of the solar convective zone using 
the helioseismology method. Gough & Sekii (1993) reported 
that they cannot find any definite evidence for the existence 
of overshooting in the Sun; while others gave an upper limit 
of 0.05-0.25Hp (Roxburgh & Vorontsov 1994, Monteiro et 
al. 1994, Christensen-Dalsgaard et al. 1995, Basu & Antia 
1991, Basu 1997). Such results are understandable. Indeed, 
what the technique of helioseismology tests is the (adiabatic) 
sound speed in the Sun, while the sound speed is determined 
by the P-T structures. In the overshooting zone, however, 
convective flux is negligible. That is why the overshooting 
below the bottom of the solar convective zone has not been 
detected by helioseismic diagnosis. From fig. (2] it is clear 
that the relative difference in sound speed between the local 
and non-local solar models is less than 1%. Such tiny differe- 
ces were indeed detected by the inversion of adiabatic sound 
speed in helioseismology (Basu 1997) , however the observed 
abrupt increase of the adiabatic sound speed at the bottom 
of solar convective zone was not correctly attributed to the 
non-local overshooting. In fact, convective flux changes its 
sign when crossing the boundary of convective zone, becom- 
ing negative in the overshooting zone (Lc < 0); in there the 
radiative flux Lr will be even larger than the total flux of the 
Sun (1/0). In the overshooting zone, the temperature gradi- 
ent will overtake the radiative counterpart V > Vrad (see 
fig. As a result, the temperature at the bottom of solar 
convective zone will rise up, just as what has been detected 
by helioseismology technique. 
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Figure 4. The super-adiabatic temperature gradient V — Vadi 
super-radiative temperature gradient V — a.nd the fractional 

convective flux Lc/L versus logP in the lower convective and 
overshooting zones for a non-local convection model of the Sun. 




Figure 6. The auto- and cross-correlations of turbulent veloc- 
ity and temperature x, Z^/^ and V versus depth (logP) for a 
non-local convection model of the Sun. The vertical dashed lines 
indicate the upper and lower boundaries of the convective zone. 
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Figure 5. A sketch of the lower convective and overshooting 
zones in the usual phenomenological non-local mixing length the- 
ory (Monteiro et al. 2000). 



Gough & Sekii (1993) measured the extension of over- 
stiooting at tlie bottom of solar convective zone following a 
picture of the overshooting zone made by the phenomeno- 
logical non-local mixing length theory, which is illustrated 
in fig. [5l In the unstable zone, the temperature gradient is 
slightly higher than the adiabatic one, while being slightly 
lower in the overshooting zone. After a distance of an over- 



shooting length dov , the temperature gradient switches sud- 
denly to the radiative from adiabatic, making a disconti- 
nuity in temperature gradient at the bottom of the over- 
shooting zone. The jump size of the temperature gradient 
is proportional to the overshooting distance dov Such a dis- 
continuity in temperature gradient is exactly what Gough 
& Sekii used to detect the overshooting distance dov It is 
the implicit assumption of full (either positive or negative) 
correlation between turbulent velocity and temperature fluc- 
tuations that makes the misunderstandings of overshooting 
zone in the non-local mixing length theory (Xiong 1985, 
Petrovay & Marik 1995). In reality, however, the turbulent 
velocity-temperature correlation decreases very quickly and 
approaches zero near the lower boundary of the solar con- 
vective zone where convective energy transfer is very effi- 
cient, as demonstrated in fig. [1] Therefore, there is no sim- 
ilarity in the structure overshooting zone between the phe- 
nomenological non-local mixing length theory and the dy- 
namic theory of non-local convection. In our view, the tem- 
perature gradient has already been smaller than the adi- 
abatic one (V < Vad < Vrad) before reaching the lower 
boundary of convective zone. The convective flux becomes 
negative passing through the boundary, therefore the tem- 
perature gradient V is smaller than the adiabatic temper- 
ature gradient V^d and higher than the radiative one Vrad 
(Vad < V < Vrad). The temperature gradient changes con- 
tinuously instead of abruptly from Vad to Vrad- The struc- 
ture of the overshooting zone in our dynamic theory of non- 
local convection is shown in fig. ID In the overshooting zone 
under the convective zone, there is a narrow (~ 0.25Hp) 
and weakly super-radiative region. Actually, the overshoot- 
ing zone is nearly radiative rather than nearly adiabatic. 
Therefore, it is not a surprise why Gough & Sekii (1993) 
could not find any firm evidence for the existence of the 
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overshooting under the bottom of solar convective zone. We 
can further justify that the methods based on stellar ther- 
mal (P-T) structure will all under-estimate the true over- 
shooting distance. In the lower overshooting, the turbulent 
velocity-temperature correlation is very small, and convec- 
tive energy transfer in there is negligible. The overshoot- 
ing distance detected by convective energy flux will be far 
smaller than that is represented by the turbulent velocity 
and temperature fields. It follows from flg. [6] that, in the 
overshooting zones either at the surface or bottom of the 
solar convective zone, the overshooting distances of turbu- 
lent velocity and temperature are both very extended, their 
e-folding lengths of overshooting are given in table [T] Our 
theoretical e-folding length agrees fairly well with those de- 
rived from observations of the solar granular fleld (Keil & 
Danfield 1978, Nesis & mattig 1989, Komm, Mattig & Nesis 
1991). 

The overshooting distance in terms of stellar evolution 
is the extension of the non-local convective mixing of chem- 
ical elements. Obviously, it is neither that of convective en- 
ergy transfer nor those of turbulent velocity and tempera- 
ture fluctuations. Calculations of massive star evolution un- 
der our complete non-local theory of convection shown that 
the non-local convective mixing overshoots a very extended 
distance (Xiong 1986). 

Although the overshooting at the bottom of the solar 
convective zone cannot be observed directly, we fortunately 
have another excellent indicator, which is the Lithium abun- 
dance of the Sun and solar type stars, for the extension of 
overshooting. No matter how disputed the mechanism of 
Lithium depletions in the atmospheres of solar type stars 
is, it can provide, at least, an upper limit for the extension 
of overshooting zone in these stars. As it is well known, '^Li 
gets burnt due to reaction ^Li{P, a)'^He at a temperature of 
2.5 X lO^K. The depth of the solar convective zone, as given 
by helioseismology, is Tc/Rq « 0.713, and the temperature 
at the bottom of the zone is Tc f» 2.26 X 10®K (Basu & 
Antia 1997, Christensen-Delsgaard et al. 1991). This tem- 
perature is not high enough to burn Lithium, without over- 
shooting (bringing Lithium deeper to higher temperatures) 
there will be no depletion of Lithium in the Sun. It is the 
overshooting that brings Lithium to the burning region at a 
higher temperature, and causes depletion. Figs. [7^-{7}l show 
the Lithium abundance depletions due to overshooting for 
M=0.90,0.95,1.00 and l.OSM© stellar models as functions 
of depth (logP), when the effect of evolution is not consid- 
ered. The dashed lines in the plots indicate the boundary of 
convective zone. As from fig. [T];, there is a gradually accel- 



erating reduction of Lithium abundance in the overshoot- 
ing zone. In the upper part of the overshooting zone, the 
mixing caused by non-local convection is very efficient, the 
abundance keeps the same as in the convectively unstable 
zone for about 0.4 Hp in length downwards. It is then fol- 
lowed by a partially mixed region where Lithium abundance 
is reduced quicker toward the center and vanishes suddenly, 
such a partial mixing region is about 0.5Hp in depth. Even 
deeper is the non-mixing zone. If the overshooting distance 
is taken as the e-folding length of the abundance from the 
bottom of the convective zone, it reads about 0.83Hp; other- 
wise if we count the deeppest bottom of the mixing process, 
it reads about 0.9Hp for the solar model. 

The equation for the conservation of Lithium abun- 
dance can be written as, 

where C is the Lithium abundance by mass, q is burn- 
ing rate of Lithium, and U is the correlation of the radial 
component of turbulent velocity w'^ and the turbulent fiuc- 
tuations of Lithium abundance defined as, 

U = ^4^, (36) 

therefore ATvr^pU is the total flux of convective mixing 
of Lithium passing through the sphere of radius r, and the 
first term on the right hand side of eq. (|35p is the rate of 
variation of Lithium abundance due to non-local convective 
mixing. It follows from fig.[6]that, in the overshooting zone, 
turbulent velocity decreases nearly exponentially. The non- 
mixed region does not mean there is no mixing at all. In our 
picture, mixing is always there, the only difference is quan- 
tity of mixing. When the non-local mixing timescale T,nix be- 
comes much longer then the nuclear timescale r„„c of deple- 
tion, i.e. the non-local mixing cannot feed fresh Lithium into 
the burning zone. Lithium abundance vanished abruptly. 
Figs [8K~d give these two terms on the right hand side of 
eq. (|35|l and their sum (the depletion rate of Lithium) versus 
depth near the lower boundary of convective zones for 0.90, 
0.95, 1.00 and 1.O5M0 main sequence stellar models respec- 
tively. The unit for these quantities are all yr~^ . The verti- 
cal dashed line locates the lower boundaries in these models. 
As clearly shown in fig. [8l the depletion of Lithium starts at 
r ~ 2.5 X 10^ K, and it goes up very quickly as temperature 
increases (being proportional to the 21st power of temper- 
ature). In the whole envelopes of these models, U < 0, this 
means that convection keeps feeding Lithium from outer into 
inner layers in order to supply the depletion at the burning 
zone. In the convectively unstable zone (where q is extremely 
small), convective mixing is very efficient, d (^ATvr^pU^ /dMr 
is nearly constant (see fig.|8]). Towards the deep interior, the 
nuclear burning rate q goes up abruptly. In the upper part of 
the overshooting zone, mixing due to convective overshoot- 
ing is still efficient enough to compensate the depletion due 
to nuclear burning, therefore the Lithium abundance profile 
is still horizontal (see figs. [7]and|8]). At the lower part of 
the overshooting zone, however, the convective overshooting 
mixing (the dotted lines in fig. [Sj can no longer support the 
balance between mixing and depletion (long dashed lines in 
fig. [HI . —^d{A-KpU)/dMr — q decrease abruptly towards the 
center, corresponding to Lithium abundance dropping off 
in fig. [T] This means the boundary of overshooting zone is 
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Figure 7. The Lithium abundance versus depth (log T) and age (labeled on the curves) for main sequence stars. The dashed vertical line 
indicates the location of the lower convective boundary. Evolution is not considered in the calculations, a). M=0.90Mq; b). M=0.95Mq, 
c). M=1.0Mq and d). M=1.05Mq. 



reached at this place. For the M=O.9OM0 star, convective 
zone is very deep with its bottom already at the burning re- 
gion of Lithium, this may leads to a shallower overshooting 
zone. Toward higher masses, convection becomes shallower, 
and the bottom of the zone goes farther away from the re- 
gion of burning (see figs.[8};-d), the overshooting zone is then 
becoming more extended. Fig. [5] presents the overshooting 
distance dov measured by Lithium depletion as a function 
of stellar mass, from which it is clear that the overshooting 
distance dov increases as stellar mass increases, going from 
0.26Hp for M=O.8OM0 to 1.65Hp for M=1.O75M0. The 
overshooting distance defined by dropping off of Lithium 
abundance by a factor of e is shown in the 8th column of ta- 



ble [T] while that defined by the distance from the boundary 
of convective zone to where Lithium becomes zero is given in 
the last column of table [1] It is clearly from figs. [7^-d that 
Lithium abundance vanishes very quickly after the e-folding 
depletion, the distance between them is less than O.lHp. 
Completely different from the Lithium abundance profile in 
fig- El turbulent velocity (x), temperature fluctuation (Z^^^) 
and velocity-temperature correlation (V) decrease exponen- 
tially with depth in the overshooting zone. The e-folding 
distances determined from the curves are given in the 2nd- 
4th colums (for the upper part of overshooting zone), and 
4th-7th columns (the lower part of overshooting zone). The 
e-folding distances given by turbulent velocity and temper- 
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Figure 8. U/C {sec burning rate of Lithium ^Li (yr ^) q, change rate of convectivc diffusion for Lithium -q = 

~'b 8M (^ff'^pC^) {yr~^) and r] — q versus depth. The dashed line indicates the lower boundary of the convectivc zone. a). M=0.90Mq; 
b). M=O.95M0, c). M=1.OM0 and d). M=1.O5M0. 



ature fields are very close to the analytic asymptote in our 
theory (Xiong 1989b). They hardly change with stellar mass, 
and are rather different from the e-folding distances defined 
by Lithium depletion. The upper and lower overshooting 
zones are a bit different in terms of these e-folding distances, 
this is due to the fact that, in the overshooting zone above 
the convective zone, gas density is low so that the effective 
Peclet number Pe <^ 1, and convective energy transfer is 
inefficient; while on the contrary, in the overshooting zone 
attach to the bottom of the convective zone, Pe ^ 1, and 
convection is highly efficient. These arguments mark the dis- 
tinct properties of the velocity-temperature correlation in 
the upper and lower overshooting zones: in the surface over- 



shooting zone, RvT ~ —1, while in the bottom overshooting 
zone, RvT ~ —0.0. When plotting fig.[T] Rvt has been mag- 
nified in order to show the details. In fact, we should have 
— 10'' < Rvt < in the lower overshooting zone, so that it 
is completely buried in the solid line of (Lc/L). 



4 SUMMARY AND DISCUSSIONS 

Detailed discussions on the definition of the boundary of 
convective zone, and the distance of convective overshoot- 
ing are presented in this paper. The main results can be 
summarized in the following: 
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Figure 9. The overshooting distance as a function of stellcir mass, 
as derived from Lithium depletion in solar type stars. 

(i) Choosing the place where the convective flux (or 
equivalently the turbulent velocity-temperature correlation) 
changes it sign as the boundary is the most proper and con- 
venient. Where the convective flux is greater than zero is 
convcctivcly unstable zone, while it is the overshooting zone 
when the convective flux is smaller than zero. The convective 
zone defined as such not only makes the local the non-local 
convection models with the same depth of convection zone 
to have similar structures, but also to have very clear phys- 
ical meanings: convective zone is the buoyant force driving 
zone for convective motion, while the overshooting zone is 
the dissipation zone against convective motion, which can 
only be supported by non-local convective difi^usion; 

(ii) It is not quite right to talk about a general overshoot- 
ing distance for stellar convection. The distance of over- 
shooting is different for different physical quantities. The 
effect of overshooting for convective energy transfer, for in- 
stance, is not important. However, the overshooting dis- 
tances of turbulent velocity and temperature fluctuations 
arc quite extended, the e-folding lengths can reach 0.25- 
0.5Hp. The overshooting distance in term of stellar evolu- 
tion is the extent of convective mixing of chemical elements. 
From the example set by the depletion of Lithium in so- 
lar type stars, it is found that convective mixing of matter 
happens in stellar evolutionary (nuclear) timescales, and is 
very efficient. Very extended and weak overshooting can still 
induce fairly efficient mixing in a very long timescale of evo- 
lution. Therefore we anticipate a very extended overshooting 
for mixing of matters, the e-folding lengths of which is gener- 
ally larger than that of turbulent velocity and temperature. 
The one in massive stellar model, for instance, can reach 1 
pressure scale height (Xioiig 1986). 

(iii) The problem of Lithium depletion in solar type stars 
is a special case of overshooting mixing that moves Lithium 
to the burning region from surface. Under such circum- 
stance, overshooting distance depends on the location of 



convective zone relative to that of Lithium burning region. 
Obviously, such conclusion cannot be simply generalized to 
the case of the core nuclear reaction. For Lithium depletion 
in solar type stars, the convective zone becomes shallower 
for higher masses, and Lithium depletion becomes slower 
(longer timescale), and this makes the extended overshoot- 
ing tail to be efficient in mixing, and the mixing range to be- 
come more extended. As a result, the overshooting distance 
detected by Lithium depletion in solar type stars becomes 
larger for increasing stellar mass. In fact, this indicates that 
mixing of matter becomes more extended for increasing nu- 
clear burning timescales. 
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